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Abstract. We demonstrate a set A and a value of s for which the sequence of W-point discrete minimal Riesz 
.s-energy configurations on A does not have an asymptotic distribution in the weak-star sense as N tends to infinity. 



1. Introduction 

In the study of classical electrostatics the representation of a large collection of electrons on a conductor 
by a charge density function is commonplace. In many cases such a representation is both physically and 
mathematically grounded. Here we consider an interaction potential derived from electrostatics and certain 
conductor geometries as a means to establish limits on when such a representation of many point charges is 
valid. In particular we demonstrate that, for certain fractal conductor geometries and Riesz potentials, there is 
no limiting charge distribution. 

Consider a conductor A as a compact subset of R'' with Hausdorff dimension d. Let co^ - {xi,. . ., x^] denote 
the locations of electrons on A. The electrostatic energy of llin is, up to a constant, 

i=i j^i 1^' 

where s is chosen to be one. By varying s the Riesz i-kernel \x - y| ' can represent generalizations of the 
Coulomb potential that decay or are singular to varying degrees. 
The infimum of the A^-point s-energy is denoted by 

&siA,N) := inf EAcon)- 

We extend £.5 so that &s{A,Q) = &s{A, 1) :- 0. For convenience we exclude the trivial case that A has finitely 
many points. For any positive value of s the functional Es is lower semicontinuous, therefore, by the compact- 
ness of A, there is at least one configuration, which is denoted <y^, that satisfies 

SAA,N)^EAc^'^), 

Identifying a minimal configuration for even simple conductor geometries such as is a formidable task. 
Alternatively, one may study qualitative properties of as tends to infinity. A common way to do this 
recasts the problem in terms of measures. For each A^ and s construct the probability measure 



xeoj' 



that places a scaled Dirac-measure 6 at each of the points in w^. One then considers measure-theoretic properties 
of the sequence of measures {yu^'^^l^^j- 

In the case when s < d one may formulate a continuous version of this problem as follows: Let Ai(A) denote 
the (unsigned) Borel measures supported on A. Let Mi(A) c M(A) denote the Borel probability supported on 
A. If we represent a charge density by a measure fi e M{A) then the Riesz s-energy of fi is 



I I - — ■^diu{y)dfi{x), 



Los Alamos National Laboratory, an affirmative action/equal opportunity employer, is operated by Los Alamos National Security, LLC, 
for the National Nuclear Security Administration of the U.S. Department of Energy under contract DE-AC52-06NA25396. 

1 



2 



MATTHEW T. CALEF 



and the electrostatic (s - 1) energy of yu is Iii/S). We may then consider the minimization problem 

yu* = ai-g mm{I,(ij) : e Mi(A)}. 
It is well known (cf ||8]) that fj.^ exists and is unique, that 

and that for any function / that is continuous on A 

The last condition is referred to as weak-star convergence of /i' '*' to A natural interpretation is that yu' is 
the continuous equilibrium charge distribution on A and that it is a limit, as tends to infinity, of minimal 
i-energy A^-point configurations. These results have practical value as it is often substantially easier to obtain 
finite-dimensional approximations of fi^ than it is to determine the points that make up for large A^. 

In the case when s > d every non-zero Borel measure supported on A has infinite i-energy (cf. |9|) and 
there is no obvious candidate for a continuous charge distribution that is the limit of the discrete minimal energy 
configurations. For this range of s recent results in |2l[3]|2l show that when A has certain rectifiability properties, 
the A^-point minimal configurations are asymptotically uniform, i.e. {jt/*''^}^^2 converges in the weak-star sense 
to a suitably normalized c/-dimensional Hausdorff measure l-C' restricted to A. Additionally, the order of growth 
of &s{A, N) is shown to be A^i+'/'' when s > d and A^^ log A^ when s = d. 

In im a candidate for a normalized s - d energy is presented. (An alternate normalization for the case 
s = d \s presented by Gustafsson and Putinar in |6|. Putinar uses this normalization to study inverse moments 
problems in [llj.) The work in |5, 4] was motivated by the hope that one could formulate, at least in the case 
s - d,?L normalized energy that could be used to obtain results for the limit as A^ tends to infinity of the discrete 
minimal tZ-energy configurations. This effort was able to show that, for certain rectifiable or fractal sets, the 
normalized energy is uniquely minimized 'H'^ restricted to A and normalized to have mass 1 - i.e. the uniform 
measure - and that converges in the weak-star sense to this uniform measure as s approaches d from below. 
However, it is not yet clear what can be inferred about the limit as A^ grows of the discrete minimal cZ-energy 
configurations. 

Relatedly, the results presented in |7| rely heavily on characteristics - e.g. the local structure of the configu- 
rations - that are lost in the weak-star limit. Further, results from |2| show that for s sufficiently large and A 
in a class of self-similar fractals these local characteristics can cause the discrete minimal energy to oscillate to 
leading order as A^ tends to infinity. These results suggest that there are limits to what can be obtained from any 
continuous minimization problem, and this paper affirms this. By using the oscillations in energy demonstrated 
in O, we demonstrate a value of s and a set A where any sequence of minimal A^-point i-energy configurations 
does not converge in the weak-star sense. For such a set there can be no valid continuous representation of the 
limiting distribution as the number of points grows to infinity. 

1.1. Main Results. We consider, for some compact K c W of Hausdorff dimension d, the functions (cf. fT^) 

,. . , &s{K,N) _ &s{K,N) 
p (A) = limmf ; — -1— and j(/r) = hmsup ■ — n— ■ 

\f g JK) - g^j(K) we denote the common value by gs ,i(K). With this we prove the following theorem. 

Theorem 1.1. Let A 'ZM.'' be the disjoint union of two compact sets A \ and A2 of Hausdorff dimension d > 
satisfying the following for some s > d: 

1. diam(Ai) and diam(A2) are both less than d{A\,A2) '■- 'vi\f{\x -y\: x e Ai, y e A2), 

2. Q<g^J,Ax)<g, j{A2) < 00, 



3. gsAAi) exists and is positive and finite. 

>N=2 



Then the sequence of measures [fi"''^]'^- C A\\{A) cannot converge in the weak-star topology on M\(A). 
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The central idea behind the proof of Theorem ll.ll is that, if &s(Ai,N) oscillates to leading order, i.e. A^'+'/'', 
then the ratio of the number of points in wj^ n Ai to the number of points in n A2 cannot be constant. This 
is sufficient to show that if one chooses a Urysohn function that is 1 on A 1 and on A2, then the limit 

lim f (f>du''^ 

N^co J 

cannot exist, which is enough to show that {//'•'^l^^j cannot have a weak-star limit. 

The rest of this paper is organized as follows: Section |2] describes the class of sets that we consider and 
provides an example in this class based on the results in Q. A particular weak-star cluster point as tends to 
infinity for the discrete minimal A^-point energy is demonstrated. Section [3] proves a lemma regarding the rate 
of change of the ratio of the minimal A^-point energy divided by its leading order in A^. Section|4]uses the results 
from the previous sections to show that, for the ranges of s and sets under consideration, there is no weak-star 
limit to the sequence {jt/'-'^l^^j. 



2. The Set A and a Weak Star Cluster Point 

We shall consider a set A that is the disjoint union of two compact set Ai and A2. We require that diam(Ai) 
and diam(A2) are both less than d{Ai,A2) = inf{|x - y\ : x e Ai,y e A2). We shall further restrict ourselves to 
such sets satisfying < ^ j*-"^'-* ^ Ssdi'^i) < °° ^^'^ < gs,d{A2) < 00. 

2.1. An example of such a set. In [7 | Hardin and SafF show that if A2 is a compact subset of a li-dimensional 
C'-manifold embedded in R'^ satisfying < '7Y''(A2) < 00, then < ^ (A2) - g^j(A2) < 00. In Proposition 
2.6 of im Borodachov, Hardin and Saff show that if Ai belongs to a certain class of self-similar (/-dimensional 
fractals, then < g ^(Ai) < g^j(Ai) < 00. The set Ai belongs to this class if 



(1) Ai=Q^,(Ai), 

1=1 



where each (p, : W —> W is a similitude with scaling L and where ^,(^1) <^ 'Pj(Ai) - for all / j. It is 
a consequence of results by Moran presented in |10| that there is a unique compact set Ai that satisfies the 
condition in Equation ([T]|, that the Hausdorff dimension of Ai denoted by d is the solution of the equation 
1 = KL'' and that < l-C'iAx) < 00. Cantor sets are a prominent example from this class. 
For a concrete example consider the set that is the union of 

A2 = [(xuxi) e : xi e [3,4] andx2 = 0} 

and the set Ai satisfying 

4 

Ai = y^,.(Ai), 

1=1 

where 

= ^ .^2W = ^ + (|oj, ^3(x) = ^+(o,^J and ^4(x) = ^+(^,^j. 
The dimension of both A\ and A2 is 1, and both have positive and finite "H' measure. 

2.2. A weak-star cluster point of {//^•'^l^^j M\{A). Here we shall identify a cluster point in the weak-stai- 
topology on A1i(A) of the sequence of measures {yu'-'^j^^j- This cluster point is one in which the .s-energy of 
the sequence of configurations {w^ n A i )^^2 achieves g ^(A 1 ). We proceed by drawing on ideas and techniques 
developed by Hardin and Saff in |7 1. 

An upper bound for ^^^J^/d obtained as follows: Choose natural numbers Mi and M2 - N -My and 

consider a configuration of points cbjsi such that n Ai - M\ and #0)^ n A2 = M2, and where EsifiiN n Ai) = 
&s{A\,Mi) and £^{(1)^ n A2) - fi.5(A2, M2). Here the # denotes the number of elements in the set following it. 
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Bounding from above the interaction energy of points on Ai and points on A2 by t/(Ai, A2) and dividing by 



Let jA^^I _j c N be an increasing sequence so that 



&s{AuN^,) 



TT77JT~ " S ,iA\). 

n— >ot) —s,d 



For Q- € (0, 1) let {N"}'^^i c N be the sequence defined by A^" - [^A^^J- Applying the bound given in Inequal- 
ity (O where is chosen to be A^", Mi is chosen to be A^^' and M2 is chosen to be A^" - A^ gives 



< 



Nl \'^'" £, (Ai,A^',) (N-^ - Nl]'"" Ss {A2,N^ - A^',) 



j^ll+s/d 



i+s/d 



For every e > we may find an A^o = Nois) sufficiently high so that for every A^" > A^o 

(3) < a^-%/A.) + (1 - ay-'^gUA2) ^ s. 

The unique value of a that minimizes the right hand side of Inequality ^ is 

gUAi)'^' 



a = 



g (Aiy/-^ + g,j(A2)'"'' 



-s,d 



Define {A^f by A^f = [^A^^J. 



To obtain a lower bound for ^'i'^;^ begin by defining the function A^i : N ^ N by 

Ni(N) = min nAi). 

Because an A^-point .s-energy-minimizing configuration may not be unique, we take a minimum over all 
A^-point .s-energy-minimizing configurations to determine the value Ni{N). We define N2{N) to be A^ - Ni{N). 
Note that A^2(A^) > min n A2). 

Given A' and a minimal A^-point configuration wj^ such that A^i (A^) = #dLt^ n A 1 , we may bound from below 
£j(A, A') = Es{cjI/) by discarding the interaction energy between points in n A\ and points in n A2 and 
then further replacing the points in n A\ with a minimal A^i (A^)-point configuration in Ai and replacing the 
points in n A2 with a minimal A^2(A^)-point configuration in A2. The bound is then 

fi.(A,A^) ^ iNmV^'''' SAAuNdN)) ^ ( N2{nW ^'" &AA2,N2(N)) 

We would like to improve the bound above by employing the asymptotic properties of ' ' ^^^^'^^J^ and 

£j(Ai,A^i(A^)) ' 
A^ • ^^^^ require ensuring that both Ni(N) and N2{N) tend to infinity as A^ tends to infinity. 

This will be accomplished by appUcations of Lemma ITTl which employs ideas presented by Bjorck in 11] . 

Lemma 2.1. Let B\, B2 c W be compact and of dimension d > 0. Further suppose diam(B2) < d{B\,B2). Let 
s > d and define 

Ni(N)^ min n Bi), 

where the minimum is taken over all N-point minimal s-energy configurations within Bi U £2- Then 

lim inf A^i (A^) = 00. 



A SEQUENCE OF DISCRETE MINIMAL ENERGY CONFIGURATIONS THAT DOES NOT CONVERGE IN THE WEAK-STAR TOPOLOGY 5 



Proof. For sake of contradiction assume that 

lim inf A^i (A^) = L < oo, 

then there is an increasing sequence {A^,)J^j c N so that N\{Ni) - L for all / e N. The contradiction shall be 
that, for large enough A^, a configuration that only places L points on B\ cannot have minimal energy. 
We first show that the quantity 

R - inf max d{y, K) 
K<zBi y"^' 

is positive. We may choose a sequence of L-point configurations in Bi, {K„]'^^^, so that 

lim v!\axd{y,K,i) - R, 

and because the L-fold product of Bi with itself is compact, we may choose {^nlj^j to be convergent to some 
K*. Continuity of the function B'^ 3 K ^ maXyeB, d{y, K) eR allows us to conclude 

m£Lxd{y,K*) = R. 

yeBi 

The set Bi is of positive Hausdorff^ dimension so it is infinite, which is sufficient to conclude R > 0. 

If NiiNi) = L for all /, then for every / there is a point r, 6 Z?i that is separated from by at least R for any 
minimal A^,-point configuration c Bi U B2. We may bound the potential energy at r, due to the point in 
from above by 

f^^'"') = ^ + W - L)d{BuB2r, 

where the first term is an upper bound for the contribution of the L points of n Bi and the second term is an 
upper bound for the A^, - L points of n B2- 

Alternatively, given any point Xj e oj-J^ n B2 we may bound from below its point energy due to the other 
points by 

^Mj^ = ^ W -L-l) diam(B2)-^ 

In this lower bound we have excluded the contribution to the point energy at xj from the points in n Z?i, and 
bounded from below the contribution to the point energy at xj from any point in n B2\{xj] by diam(B2) 

Because B2) > diam(B2)^ we may find A^,- sufficiently high so that t/(r,) < U jixj), but then the energy 
of configuration could be reduced by moving the point from Xj to rj and this contradicts the assumption 
that <jj^ has minimal energy. □ 



We may apply Lemma [ZTI to the sets under consideration by identifying either Ai or A2 as Bj. This is 
sufficient to show that 

lim inf A^i (A^) - liminf A^2(A^) - ^■ 

N—^co N—>co 

Along the the sequence {Af"')^[, we may apply the bound given in Inequality (O and, because Ni{N" ) and 
A^i (A^^ ) grow to infinity, we may apply Inequality (|4|i giving, for A^" sufficiently high 

£,(Ai,A^,(A?,f)) £.(A2,A^2(A^,f)) 
Here we have bounded from below — — — by? (Ai)-e/2 and boundedfrom below — — — 

N^m') ' N2m') ' 

hygsAA2)-sl2. 

Because e is arbitrary and because a* is the unique minimizer of the left hand side of the above upper bound 
we may conclude that 

hm : — - a . 
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If we let if/ 6 All (A) be any weak-star cluster point of [j/'^" Y^^^, and if we choose a continuous function 
; A -> R such that <^(jc) = 1 for all ;c e Ai and (p{y) - for all y e A2, then 



<pd>f/ - ip{A\) - a*. 
This is sufficient to prove Lemma lZ2l 

Lemma 2.2. Let A be the disjoint union of two compact sets A \ and A2 that meet the following conditions: 
1. both Ai and A2 are of Hausdorjf dimension d. 



2. diam(Ai) ami diam(A2) are both less than c/(Ai, A2), 

3. for some s > d, < g JAi) < ^.5,/(Ai) < 00 and 

4. for the same s, gs,d(A2) exists and < gs,d(A2) < ^■ 



Then there is a weak-star cluster-point i// 6 Mi(A) of the sequence of measures {yu*''^}^^, ^^^^ 



mi) 



g (Ai)-\- g,AA2) 

—s.d 



3.RateopChanoe™^4^pok^cR.Compact 

In this section we consider a compact set K izW so that < f , J^K) < 00, and bound from below the rate of 
change of the function G : N ^ R defined by 

This is the content of Lemma ITTI In Corollary [32] we use Lemma [TTI to show that, if for some A^o, G(No) is 
close to g^j(K), then for some A^' that is larger than A^o and for which the ratio Nq/N' is close to 1, we will have 
that G{N') is also close to g^ d(^)- 

Lemma 3.1. Let K <z W be compact and s > d — dim/T > 0. Let N > 2 and N' > I be natural numbers and 
let K = j^. Then 

G((l + k)N) > G{N) -(1 + ^)'<G(N). 

Proof. From our assumptions 

I , \I+.s/rf 

- (ttt) G{N) 
^ (i-(i + 3)^)g(a^) 

= G(Ar)-(l + |)^G(A^), 

/ . \l+s/£/ 

where the last inequality follows from a first-order expansion of h{K) '--XytA about zero and the convexity 
of h. □ 

Corollary 3.2. Let K (zW be compact and s, d > so that < gsjiK) < 00. Let {A/„}^j cN be an increasing 
sequence so that 

&AK,M„) _ 

For every e > there is a 6 > and an Mq < 00 so that if for some M' E N and M e {M„}J^j, 

M 



Mq < M < M' < 



1-5' 
then 

&AK,M') 
^,1,,;, >8s,diK)-s. 
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Proof. Let e > be arbitrary. Choose Mq so that for any M e {M„)JJ1, greater than Mq 
Choose /<•() sufficiently small so that for all k e (0, kq), 

&, (K, M) s 

(6) (1 + s/d)K sun — - — -— < -. 

M>M„ M^+^/d 2 

Combining Lemma [3T| and Inequalities (|5]l and (|6]) gives for any M e N and k e (0, a^q) where (1 + k)M € N 

£j iK,il+ k)m) (k, m) (k, m) 

(7) — r-^ > — --(l + s/d)K^ ->g,JK)-e. 



Choose 6 . If M' and M are such that Mq < M < M' < j^j, and if M' = (1 + k)M, then k < and the 

bound in Inequality ^ ensures 



&s{K,M') ^_ 

>gsAK)-e. 

□ 

4. NoN-CoNVERGENCE OF {/i^ '^f^^j ™ THE WeAK-StAR ToPOLOGY 

In this section we prove Theorem ll.il - that the sequence of measures [fi^-'^]'^^^ c M\{A) cannot converge 
in the weak-star topology on M\{A). If the sequence did converge, it must converge to the measure ijj identified 
in Lemma l272l this will be shown to lead to a contradiction. 



Proof of Theorem For sake of contradiction, assume that the sequence [fi^''^}'^^^ c M\(A) converged in the 
weak-star topology on M\{A). Since this topology separates elements of M\{A) and since, by Lemma IZTl a 
subsequence of {;u''^)^^2 converges to i/r, the assumption would require that ^"'^ converges to i/r in the weak-star 
sense as — > oo implying 

(8) hm — - — - a - 



N g iAi)+gs,d{A2) 



Let |-^n| c N be an increasing sequence so that 



Our intention is to find A^'" E |^!| an A^' e N so that we may apply CoroUarv 13.21 where M is 

identified with N'" and M' is identified with N](N'). 

We proceed as follows: Let e > 0. Let 6 and Mq be as provided by Corollarv l3 .21 applied to the case K = Ai 

f— 1 1 °° , 2y 
and {M,,)"! J — lN„ \ . One may verify that it is possible to choose y e (0, 1) such that, for y' , the 

[ )n=i 1 - 7 

following (motivated by inequalities arising later in the proof) hold: 

I \ . . 1 + 

We choose such a y. Define the non-decreasing sequence {Nn}'^^^ by the equation N„ = (1 + 3y)Nl 
From Equation (O and from Lemma l2!2l we may further increase Mq so that the following all hold 

NdN) 



7' + (l+7') <<5 and 7-^ > 1- 

1+37/ 1+7' 



(9) sup 

N>Mo 



N 



< min{7a*,e}. 



g 
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(10) 

(11) 

Inequality ^ implies 
(12) 



&AA2,N2(N)) 



and 



sup 

N>Mo 



N2(N) 



l+s/d 



- gsAAi): 



< E. 



1 + -y I < L(l + 3r)A^J for all > Mq. 



sup 

N>Mo 



- 1 



Choose A^' e N and A^" € so that both are greater than Mo, and so that 

N" 



(13) 

This allows the following bound 

(14) 



a*N' 



<y. 



N" 



Ni(N') 



1 



< |1 

< 

< 



N" . N" _ N" 
a'N' a'N' N,(N')\ 

I N" N" 



N" 
a'N' 



a'N' N,(N')\ 

y + (jll.)\l _ JilALl 

'^\a'N')\^ Nt(N')\ 



r + (1 +r)T^ 



The last inequality above follows from Inequalities (fT2l) and ( fTSl ). Now let A^'" denote the element of |^h| 
so that A^" = L(l + iy)N"'i. Then 

N'" 



Ni(N') 



- 1^ ' 
< 



N" 1 N" 


N'" 


Ni(N') N,(N') 


N,(N')\ 


N" III N" 


N'" 


N,(N')\ ' l^iCA") 


N,(N') 


y'+ ^" ll 

' N,(N') 1 


N'" 1 
N" 1 



< 
< 



N'" 



r' +(!+/) |1- -i+3rJ^"' 

r' + (i+r')(i-TT3^). 



Here the second to last inequality follows from Inequality (fl4l l. Inequality (fl4] l also implies 

A^" A^" 1(1+ 3r)A^"' I (l + f t) 

1 - — > -y, hence Ni(N') > —— = ^ > N'"^ ^ 

A/i(A^') ' y' + l y' + l y' + \ 

Because y was chosen to ensure that 

^ ^ > 1, and 7' + (l+y') l \<S 



y' + l \ 1 + 3^ 

we have that Ni(N') > N'" and 

1 < 6. 

Ni(N') 

With this we may employ Corollarv l3.2| where M is identified with N'" and M' is identified with N-i{N'), giving 
the bound 
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This, combined with InequaUties (HI and (fTTT i. gives 



^/i+.s/rf ~ \ N' I ' \ N' 

In light of InequaUty ^ we have 

^'j^ta? ^ - ^)y"'"(gs,AAi) -s) + (l-a*- sy^'^'igUAi) - s). 

In the preceding arguments, e was an arbitrary positive number that constrained A^' to be above some value. 
We may then construct a sequence (A^/,}^; of such A^' so that 

lim ^ii^ > a^'^^^/A,) + (1 - a )'"/V..(A2). 

Now, for each e {A^/,}^7=i ^^^^ larger than 2M(), and for any two natural numbers M,| and M^-N'^^- Ml 
both greater than Mo, Inequalities (l2]i, (fTOl t and (fTTT t give 



For an arbitrary /? e (0, 1) we may choose and M}^ so that jfr ^ /3 and — > l-y8asn— >oo. Note also 

that, because goes to infinity as n goes to infinity, the term A^^'"'^''(i(Ai, A2) goes to zero as n grows. For any 
j6 we obtain the following bound 

(15) /^''■'"g.AAi) + (1 -/Sy^"'gsAA2) > lim > a*'''"g,JAi) + (1 - a )'-'/'^g.,,(A2). 

The function f(J3) - j8'^ '^''l.5,rf(Ai) + (1 - yS)'+'/''gj,d(A2) is uniquely minimized by 

gsAAi) ^ gsAAi) _ ^. 
gs.d(A 1 ) + §^(^2) g^ JAi)+ g.UA2) 

invalidating Inequality ( fTsT l. This is sufficient to show that {//'•'^j^^j cannot converge in the weak-star topology 
on All (A). □ 
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